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Abstract This paper is concerned with the adaptive numerical 
treatment of stochastic partial differential equations. Our method of 
choice is Rothe’s method. We use the implicit Euler scheme for the 
time discretization. Consequently, in each step, an elliptic equation 
with random right-hand side has to be solved. In practice, this cannot 
be performed exactly, so that efficient numerical methods are needed. 
Well-established adaptive wavelet or finite-element schemes, which are 
guaranteed to converge with optimal order, suggest themselves. We 
investigate how the errors corresponding to the adaptive spatial dis¬ 
cretization propagate in time, and we show how in each time step the 
tolerances have to be chosen such that the resulting perturbed dis¬ 
cretization scheme realizes the same order of convergence as the one 
with exact evaluations of the elliptic subproblems. 
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1 Introduction 

This paper is concerned with the numerical treatment of stochastic evolution 
equations of the form 

du(t) = (Au(t) + f (u(t))) dt + B(u(t)) dW(t), u(0) = uq, (1) 

on the time interval [0, T] in a real and separable Hilbert space U. Here, 
A : D(A) C U —>• U is a densely defined, strictly negative definite, self- 
adjoint, linear operator such that zero belongs to the resolvent set and the in¬ 
verse A _1 is compact on U. The forcing terms / : D((—A) Q ) —>• D((—A) Q ~ a ) 
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and B : D((—A) e ) —>• £(£ 2 , D((—A) 8 are Lipschitz continuous maps for 
suitable constants g, a and /3; and finally, W = (IU(f)) tg [o,T] is a cylindrical 
Wiener process on the sequence space £2 = ^(N). In practical applications, 
evolution equations of the form (1) are abstract formulations of stochastic 
partial differential equations (SPDEs, for short): Usually A is a differential 
operator, / a linear or nonlinear forcing term and B(u(t))dW(t) describes 
additive or multiplicative noise. They are models, e.g., for reaction diffusion 
processes corrupted by noise, which are frequently used for the mathemati¬ 
cal description of biological, chemical and physical processes. Details on the 
equation, the operators A, the forcing terms / and B and the initial condi¬ 
tion uq are given in Section 2. Usually, the exact solution of (1) cannot be 
computed explicitly, so that numerical schemes for the constructive approxi¬ 
mation of the solutions are needed. For stochastic parabolic equations, there 
are two principally different approaches: the vertical method of lines and the 
horizontal method of lines. The former starts with an approximation first 
in space and then in time. We refer to [23-25, 27] for detailed information. 
The latter starts with a discretization first in time and then in space; it is 
also known as Rothe’s method. In the stochastic setting, it has been studied, 
e.g., in [6,22]. These references are indicative and by no means complete. 

Very often, the vertical method of lines is preferred since, at first sight, it 
seems to be a little bit simpler. Indeed, after the disretization in space is per¬ 
formed, just an ordinary finite dimensional stochastic differential equation 
(SDE, for short) in time direction has to be solved, and there exists a huge 
amount of approaches for the numerical treatment of SDEs. However, there 
are also certain drawbacks. In many applications, the utilization of adaptive 
strategies to increase efficiency is almost unavoidable. In the context of the 
vertical method of lines, the combination with spatial adaptivity is at least 
not straightforward. In contrast, for the horizontal method of lines, the fol¬ 
lowing natural approach suggests itself. Using Rothe’s method, the SPDE 
can be interpreted as an abstract Cauchy problem, i.e., as a stochastic differ¬ 
ential equation in some suitable function spaces. Then, in time direction we 
might use an SDE-solver with step size control. This solver must be based 
on an implicit discretization scheme since the equation under consideration 
is usually stiff. Consequently, in each time step, an elliptic equation with 
random right-hand side has to be solved. To this end, as a second level of 
adaptivity, adaptive numerical schemes that are well-established for deter¬ 
ministic equations, can be used. We refer to [9,10,15] for suitable wavelet 
methods, and to [13,5, IT 20,28,36,37] for the finite element case. As before, 
these lists are not complete. 

Although this combination with adaptive strategies is natural, the math¬ 
ematical analysis of the resulting schemes seems to be still in its infancy. In 
the stochastic setting, Rothe’s method with exact evaluation of the elliptic 
subproblems, has been considered, e.g., in [6,22], and explicit convergence 
rates have been established, e.g., in [12,13,26,35]. First results concern- 
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ing the combination with adaptive space discretization methods based on 
wavelets have been shown in [31]. 

Even for the deterministic case, not many results concerning a rigorous 
convergence and complexity analysis of the overall scheme seem to be avail¬ 
able. To our best knowledge, the most far reaching achievements have been 
obtained in [7]. In this paper, it has been clarified how the tolerances for the 
elliptic subproblems in each time step have to be tuned so that the overall 
(perturbed) discretization scheme realizes the same order of convergence (in 
time direction) as the unperturbed one. Moreover, based on concepts from 
approximation theory and function space theory, respectively, a complexity 
analysis of the overall scheme has been derived. It is the aim of this paper to 
generalize the analysis presented in [7] to SPDEs of the form (1). We mainly 
consider the case of the implicit Euler scheme, and we concentrate on the 
convergence analysis. To our best knowledge, no result in this direction has 
been reported yet. Complexity estimates are beyond the scope of this work 
and will be presented in a forthcoming paper. 

For reader’s convenience, let us briefly recall the basic approach of [7] for 
the deterministic case, confined to the implicit Euler scheme. As a typical 
example, let us consider the deterministic heat equation 

u'(t) = A u(t) + /(£, u(t)) on O, t € (0, T ], 
u = 0 on dO , t € (0, T], 

u( 0) = uq on O, 

where O C M rf , d > 1, denotes a bounded Lipschitz domain. We discretize 
this equation by means of a linearly implicit Euler scheme with uniform time 
steps. Let K G N be the number of subdivisions of the time interval [0, T], 
The step size will be denoted by r := T/K, and the fc-th point in time is 
denoted by t k := rk , k € {0,... , K}. The linearly implicit Euler scheme, 
starting at uq, is given by 


V j k +1 ^k 

T 


Au k+1 + f{t k ,u k ), 


i.e., 


(■I - rA)u k+1 = u k + rf(t k , u k ), 


for k = 0,..., K — 1. If we assume that the elliptic problem 


( 2 ) 


L r v := (7 — tA)v = g on O, v\qo = 0, 
can be solved exactly, then one step of the scheme (2) can be written as 

^k+l ^T,kiV'k): ( 3 ) 


where 

Rr,k(w) := W + rf(t k ,w) 
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and L t is a boundedly invertible operator between suitable Hilbert spaces. 

In practice, the elliptic problems in (3) cannot be evaluated exactly. 
Instead, we employ a ‘black box’ numerical scheme, which for any prescribed 
tolerance e > 0 yields an approximation [u] e of v := L~ 1 R rk (w), where w is 
an element of a suitable Hilbert space, i.e., 

||w - Hell < 

for a proper norm ||-||. What we have in mind are applications of adaptive 
wavelet solvers, which are guaranteed to converge with optimal order, as de¬ 
veloped, e.g., in [9], combined with efficient evaluations of the nonlinearities 
/ as they can be found, e.g., in [11,16,30]. In [7] we have investigated how 
the error propagates within the linearly implicit Euler scheme and how the 
tolerances e k in each time step have to be chosen, such that we obtain the 
same order of convergence as in the case of exact evaluation of the ellip¬ 
tic problems. We have shown that the tolerances depend on the Lipschitz 
constants of the operators 

E r .j,k Hr Rr,k— l) ° (-^r Rr,k— 2 ) O • • • O (_L r 72 t j), 
with 1 < j < k < K , K € N, via 


fc-i 

u{t k ) - Uk\\ < | \u(t k ) - u k II + Y. CyJ+i.fc £j, 

3 =0 


where u k is the solution to the inexactly evaluated Euler scheme at time t k . 

Now let us come back to SPDEs of the form (1). Once again, for the 
(adaptive) numerical treatment of (1) we consider for K € N and r := T/K 
the linearly implicit Euler scheme 


u k +1 = (I ~ tA) 1 (u k + rf(u k ) + \[t B(u k )xk) , 

k = 0,..., K — 1, 


( 4 ) 


with 

X* := Xf := ^ {W{t« + 1) - W{t £)) , 

where t k := rk. k = 0,..., K. If we set 

Rr,kH '■= W + Tf{w) + \prB(w)x k , k = 0,..., K - 1, 

L~ l w := (/ — rd) _1 w, k = 1,..., K, 


the operators being defined between suitable Hilbert spaces 7i k and Q k , the 
scheme (4) can again be rewritten as 

u k+1 = L~ lR T}k (u k ), k = 0,... ,K - 1. (5) 
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We refer to Section 3 for a precise formulation of this scheme. 

Once again the elliptic problems in (5) cannot be evaluated exactly. 
Similar to the deterministic setting, we assume that we have at hand a 
‘black box’ numerical scheme, which for any required w approximates 

v := (I ~ tA)~ 1 (w + f(w) + B(w)xk) 

with a prescribed tolerance e > 0. What we have in mind are applications 
of some deterministic solver for elliptic equations to individual realizations, 
e.g., an optimal adaptive wavelet solver as developed in [9], combined with 
proper evaluations of the nonlinearities / and B, see, e.g., [11,16,30], and 
an adequate truncation of the noise. It is the aim of this paper to inves¬ 
tigate how the error propagates within the linearly implicit inexact Euler 
scheme for SPDEs (cf. Proposition 4.3) and how the tolerances in each 
time step have to be chosen, such that we obtain the same order of conver¬ 
gence (in time direction) for the inexact scheme as for its exact counterpart 
(cf. Theorem 4.2). 

Concerning the setting, we follow [35] and impose rather restrictive con¬ 
ditions on the different parts of Eq. (1). This allows us to focus on our main 
goal, i.e., the analysis of the error of the inexact counterpart of the Euler 
scheme (4), without spending too much time on explaining details regarding 
the underlying setting, cf. Remark 2.12. Compared with [35] we allow the 
spatial regularity of the whole setting to be ‘shifted’ in terms of the addi¬ 
tional parameter g. In concrete applications to parabolic SPDEs, this will 
lead to estimates of the discretization error in terms of the numerically im¬ 
portant energy norm, cf. Example 2.11, provided that the initial condition 
uq and the forcing terms / and B are sufficiently regular. 

A different approach has been presented in [31], where additive noise is 
considered, a splitting method is applied, and adaptivity is only used for 
the deterministic part of the equation. We remark that the use of spatially 
adaptive schemes is useful especially for stochastic equations, where singu¬ 
larities appear naturally near the boundary due to the irregular behaviour 
of the noise, cf. [8] and the references therein. 

We choose the following outline. In Section 2 we present the setting and 
some examples of equations that fit into this setting. In Section 3 we show 
how to reformulate the linearly implicit Euler scheme as an abstract Rothe 
scheme and derive convergence rates under the assumption that we can 
evaluate the subproblems (5) exactly. We drop this assumption in Section 4 
and focus on how to choose the tolerances for each subproblem, such that 
we can achieve the same order of convergence. 

2 Setting 

In this section we describe the underlying setting in detail. It coincides with 
the one in [35] (‘shifted’ by g > 0). Furthermore we define the solution 
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concept under consideration and give some examples of equations, which fit 
into this setting. 

We start with assumptions on the linear operator in Eq. (1). 
Assumption 2.1. The operator A : D(A ) C U —>• U is linear, densely de¬ 
fined, strictly negative definite and self-adjoint. Zero belongs to the resolvent 
set of A and the inverse A : t/ —>• fT is compact. There exists an a > 0 
such that (— A)~ a is a trace class operator on U. 

To simplify notation, the separable real Hilbert space U is always as¬ 
sumed to be infinite-dimensional. Under the assumption above, it follows 
that A enjoys a spectral decomposition of the form 

A y = 'y]Xj( v ,e j ) u ej, v € D(A), (6) 

je N 

where (ej)j^ is an orthonormal basis of U consisting of eigenvectors of A 
with strictly negative eigenvalues (A j)jeN such that 

0 > Ai > A 2 > ... > Xj —» — 00 , j —> 00 . (7) 

For s>0we set 

OO 

D{(-AY) := {„ € U : £ |(-A,)*(»,e,-)t,| 2 < oo}, (8) 

3 = 1 

(-A) s u := J^(-A j) s (v, ej)uej, v € D((-A) s ), (9) 

je N 

so that D((— A) s ), endowed with the norm || ■ ||_d((_a) s ) := ||(—■A) s • ||c/, is a 
Hilbert space; by construction this norm is equivalent to the graph norm of 

(-Air- 

For s < 0 we define D((—A) s ) as the completion of U with respect to the 
norm || • ||d((_a)<>), defined on U by ||u|||, (( _ A)s) := EjeN I (~ X j) S ( v > e j)u\ ■ 
Thus, D((—A) s ) can be considered as a space of formal sums 

v = v^ej, such that |(— Xj) 8 v^\ 2 < oo 

j€ N j€ N 

with coefficients € M. Generalizing (9) in the obvious way, we obtain op¬ 
erators (—A) s , s£l, which map D((—A ) r ) isometrically onto D((—A) r ~ s ) 
for all r € M. 

The trace class condition in Assumption 2.1 can now be reformulated as 
the requirement that there exists an a > 0 such that 

nM)-” = E(- A >r“ < °°- (10) 

je N 
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Note that any linear operator with a spectral decomposition as in (6) and 
eigenvalues as in (7) and (10) fulfills Assumption 2.1. Let us consider a 
prime example of such an operator. Throughout this paper, we write 7/2(0) 
for the space of quadratically Lebesgue-integrable real-valued functions on 
a Borel-measurable subset 0 of M. d . Furthermore, C(Ui\U 2 ) stands for the 
space of bounded linear operators between two Hilbert spaces U\ and 02- If 
the Hilbert spaces coincide, we simply write C{U\) instead of C{U\\U{). 

Example 2.2. Let 0 be a bounded open subset of M d , set U := 7/2(0) and 
let A := Ag be the Dirichlet-Laplacian on 0, i.e., 

Ag : D(Aq) C 7/ 2 (0) 7/2(0) 

with domain 

d rj 2 

D( Ag) = {ue 771(0) :A , /: ^_„ e 7/2(0)}, 

i=l * 

where 77 q( 0) stands for the completion in the 7/2(0)-Sobolev space of or¬ 
der one of the set Cfi°(0) of infinitely differentiable functions with com¬ 
pact support in 0. Note that this definition of the domain of the Dirichlet- 
Laplacian is consistent with the definition of 77((—Ag) s ) for s = 1 in (8), 
see, e.g., [32, Remark 1.13] for details. This linear operator fulfills Assump¬ 
tion 2.1 for all a > d/2: It is well-known that it is densely defined, self- 
adjoint, and strictly negative definite. Furthermore it possesses a compact 
inverse (Ag) -1 : 7/2(0) —)■ 7/2(0), see, e.g., [21]. Moreover, Weyl’s law states 
that 



see [4], implying that (10) holds for all a > d/2. 

Next we state the assumptions on the forcing terms / and B. 
Assumption 2.3. For certain smoothness parameters 

1 — OL 

q > 0, a < 1 and f3 < —-— (11) 

(a as in Assumption 2.1), we have 

/ : D((-A)°) -»■ D((-A)°-°), 
B:D((-Ay)^C(l 2] D((-A)^)). 

Furthermore, / and B are globally Lipschitz continuous, that is, there exist 
positive constants 0j ip and 0g ip such that for all v, w € D((—A) e ), 

I \f(v) ~ f(w)\\ D ((- A )e—) < CY P \\V - w\\ D{{ _ A)e) , 

and 

II B{v) - B(w)\\c(t 2 ;m-A)e-P)) < Cb 1P \\v ~ H|p((-A)e)- 
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Remark 2.4. (%) The parameters a and /? in Assumption 2.3 are allowed 
to be negative. 

(ii) Assumption 2.3 follows the lines of [35] (‘shifted’ by g > 0). The linear 
growth conditions (3.5) and (3.7) therein follow from the (global) Lipschitz 
continuity of the mappings / and B. 

Finally, we describe the noise and the initial condition in Eq. (1). For 
the notion of a normal filtration we refer to [34]. 

Assumption 2.5. The noise W = (W(t))^ e [o ) r] is a cylindrical Wiener 
process on l 2 with respect to a normal filtration (&t)te[o,T\- The underlying 
probability space (fl, &. P) is complete. For g as in Assumption 2.3, the 
initial condition uq in Eq. (1) satisfies 

u o £L 2 (Q,^ o ,P]D({-A) 0 )). 

In this paper we consider a mild solution concept. To this end let (e tA )t >0 
be the strongly continuous semigroup of contractions on U generated by A. 
Definition 2.6. A mild solution to Eq. (1) (in D((—A) 0 )) is a predictable 
process u : Pi x [0, T] —>• D({—A) e ) with 

sup E||u(t)|||, (( _ A)e) < 00 , (12) 

te[o,T] 

such that for every t € [0, T] the equality 

u(t) = e tA u 0 + f e^ A f(u{s)) ds + f e^ A B{u{s)) dW(s) (13) 

Jo Jo 

holds P-almost surely in D((—A) g ). 

Remark 2.7. (i) Let u : Q x [0,T] —> D({—A) Q ) be a predictable process 
fulfilling (12). Then, the first integral in (13) is meant to be a D((—A) e )- 
valued Bochner integral for P-almost every uj £ Cl] the second integral is a 
D((— A) e )-valued stochastic integral as defined, e.g., in [14,34], Both inte¬ 
grals exist due to (12) and Assumptions 2.1 and 2.3. For example, consider¬ 
ing the stochastic integral in (13), we know that it exists as an element of 
L 2 (Cl, P; D((— A) 0 )) if the integral 

/' , E|| e «-*)^B( U ( S ))||L s(felD(( _ /1 ,» )) d S (14) 

•J 0 

is finite, where £hs(-^ 2 ; D((— A) 0 )) denotes the space of Hilbert-Schmidt op¬ 
erators from t 2 to D((—A) 0 ). The integrand in (14) can be estimated from 
above by 

TM— A)- ”|K-4l)' ,+ “ /2 e< , -') A |l? :(D(( . A) . )) E|| ( -4l ) -' , B ( » (s)) |e ( , 2iD(( . /1 ,. )) . 

and we have 

IK— 
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for 2/3 + a € [0,1). For 2/3 + a < 0 we have (— Xj) 2l3+a < (—Ai) 2/3+q for all 
j € No, which yields 

||(_ j4 )|5+«/2 e «-.)^||^ D( (_ j , < „ < C . 

Moreover, by the global Lipschitz continuity of the mapping B : D((—A) s ) —>• 
C(£r,D((-A)B-P)), 

E\\(-A)-^B(u(s))\\ 2 cii2 . Da _A )e)) < C(l + sup E||u(r)||^ (( _ A)e) ). 

re [o,r] 

Thus, the stochastic integral in (13) is well-defined. 

(ii) For the case g = 0 existence and uniqueness of a mild solution to 
Eq. (1) has been stated in [35, Proposition 3.1]. The proof consists of a 
modification of the proof of Theorem 7.4 in [14]—a contraction argument in 
Loo([ 0,T];L2(H;C/)). For the general case g > 0 existence and uniqueness 
can be proved analogously, see [29, Theorem 5.1]. Alternatively, the case 
g > 0 can be traced back to the case g = 0 as described in the proof of 
Proposition 2.8 below. 

Proposition 2.8. Let Assumptions 2.1, 2.3 and 2.5 be fulfilled. Then, Eq. (1) 
has a unique (up to modifications) mild solution in D((—A) e ). 


Proof. If Assumptions 2.1, 2.3 and 2.5 are fulfilled for g = 0, Eq. (1) fits 
into the setting of [35] (the Hilbert space U is denoted by H there). By 
Proposition 3.1 therein there exists a unique mild solution u to Eq. (1). Now 
suppose that Assumptions 2.1, 2.3 and 2.5 hold for some g > 0. Set 

U := D((—A) e ), D(A) := D((-A) e+1 ) 
and consider the unbounded operator A on U given by 
A : D(A) CU — >• U, v i y Av := Av. 

Note that A fulfills Assumption 2.1 with A, D{A) and U replaced by A, 
D(A) and U, respectively. Defining the spaces D((—A) s ) analogously to the 
spaces D((—A) s ), we have D((—A) Q+S ) = D((—A) s ), s € M, so that Assump¬ 
tions 2.3 and 2.5 can be reformulated with g, D((—A) s ), D((—A) e ~ a ) and 
D((-Ay-P) replaced by g := 0, D((-i)e), D{{-Ay~ a ) and D((-A)e~P), 
respectively. Thus, the equation 

du(t) = ( Au(t ) + /(u(t))) df + B(u{t)) dW(t), u( 0) = uq, (15) 


fits into the setting of [35] (now U corresponds to the space H there), so that, 
by [35, Proposition 3.1], there exists a unique mild solution u to Eq. (15). 
Since the operators e tA € C(U) and e tA € C(U) coincide on U C U, it is 
clear that any mild solution to Eq. (15) is a mild solution to Eq. (1) and 
vice versa. □ 
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Remark 2.9. If the initial condition uq belongs to L p (Q, D((—A) Q )) 

C L 2 (H, IP; D((—A) e )) for some p > 2, the solution u even satisfies 
sup te [o iT ]E||w(f)||^((_ A)£ ,) < 00 . This is a consequence of the Burkholder- 
Davis-Gundy inequality, cf. [14, Theorem 7.4] or [35, Proposition 3.1]. Anal¬ 
ogous improvements are valid for the estimates in Propositions 3.2 and 4.3 
below. 

We finish this section with concrete examples for stochastic PDEs that 
fit into our setting. 

Example 2.10. Let O be an open and bounded subset of R d , U := 1*2(0), 
and let A = Aq be the Dirichlet-Laplacian on O as described in Example 2.2. 
We consider examples for stochastic PDEs in dimension d = 1 and d > 2. 
First, let O Cl 1 be one-dimensional and consider the problem 

du(t,x) = A x u(t,x) dt + g{u(t,x)) dt + h(u{t,x)) dW\{t, x), 

(f, x) G [0, T] x O, > 

u(t,x ) = 0, (t,x) G [0, T] x dO, 

u(0, x) = uq(x), xGO, 


where uq G L 2 (0), g : M —» R and h : R —>• R are globally Lipschitz 
continuous, and W\ = (P /r i(t))tg[o,T] is a Wiener process (with respect to a 
normal filtration on a complete probability space) whose Cameron-Martin 
space is some space of functions on O that is continuously embedded in 
Loo(O), e.g., W 1 is a Wiener process with Cameron-Martin space H s (0 ) 
for some s > 1/2. Let (ii> k ) k m be an arbitrary orthonormal basis of the 
Cameron-Martin space of W\ and define / and B as the Nemytskii type 
operators 


f(v)(x):=g(v(x)), v£L 2 (0),x£0, 

(B(v)a) (x) := h(v(x )) ^ a k i> k (x ), v G L 2 {0), a = ( a k ) km G i 2 , x G O. 
k£ N 

(17) 

Then, Eq. (1) is an abstract version of problem (16), and the mappings / 
and B are globally Lipschitz continuous (and thus linearly growing) from 
D((—A)°) = L 2 (0 ) to L 2 (0 ) and from D((—A)°) to C,{k 2 \L 2 (0)), respec¬ 
tively. For B this follows from the estimate 


\\B(vi)a - B(v 2 )a\\ L2 ( 0) = (h(v 1 ) - h(v 2 )) E 


km 


L 2 (0) 


< 


||/i(ui) - h{v 2 )\\ L2{0 )Y^a k ^ k 


km 


Loc(O) 


< C'llui — ^ 2 |U 2 (o) l|a|U 2 , 


where the last step is due to the Lipschitz property of h and the assumption 
that the Cameron-Martin space of W\ is continuously embedded in L^O). 
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It follows that Assumptions 2.1, 2.3 and 2.5 are fulfilled for 1/2 < a < 1 
(compare Example 2.2) and g = a = /3 = 0. 

Now let O C M. d be d-dimensional, d > 2, and consider the problem (16) 
where uq € L2(0), g : M —>• R is globally Lipschitz continuous, /i : R —>• M 
is constant (additive noise), and H/ = (Wi(t))t e [ 0 ,T] is a Wiener process 
whose Cameron-Martin space is some space of functions on O that is con¬ 
tinuously embedded in D((—A)~d) f or SO me /3 < 1/2 — d/4. One easily sees 
that the mappings / and B, defined as in (17), are globally Lipschitz con¬ 
tinuous (and thus linearly growing) from D((—A)°) = L2(0) to 1/2(0) and 
from D((—A)°) to C(l 2 ]D((—A)~^)), respectively. It follows that Assump¬ 
tions 2.1, 2.3 and 2.5 are fulfilled for (3 < 1/2 — d/4, d/2 < a < 1 — 2/3, 
and g = <7 = 0. Alternatively, we could assume h to be sufficiently smooth 
and replace h(u(t,x)) in problem (16) by, e.g., h( f 0 k(x,y)u(t,y) dy) with 
a sufficiently smooth kernel k : O x O — > M. 

Example 2.11. As in Examples 2.2 and 2.10, let A = Aq be the Dirichlet- 
Laplacian on an open and bounded domain O C From the numerical 
point of view, we are especially interested in stochastic PDEs of type (1) 
with g = 1/2. In this case the solution process takes values in the space 
D((-A) 1 / 2 ) = Hl(0), and, as we will see later in Proposition 3.2 and The¬ 
orem 4.2, we obtain estimates for the approximation error in terms of the 
energy norm 

IMId((—A g)V2) = (Vu, Vu)^ 2 (0) , v € Hq(O). 

The energy norm is crucial because error estimates for numerical solvers of 
elliptic problems (which we want to apply in each time step) are usually 
expressed in terms of this norm, compare [7, Section 4], where adaptive 
wavelet solvers with optimal convergence rates are considered. 

First, let O C M 1 be one-dimensional, and consider the problem (16) 
where uq € Hq( 0), g : M —>• M is globally Lipschitz continuous, h : M —>• M 
is linear or constant, and Wi = (Wi(t)) te [ 0 ^ is a Wiener process whose 
Cameron-Martin space is some space of functions on O that is continuously 
embedded in for some nonnegative (3 < 1/4, so that W\ 

takes values in a bigger Hilbert space, say, in D((— A)" 1 / 4 ). (The embed¬ 
ding D((—A) l / 2 ~P) D((— A) -1 / 4 ) is Hilbert-Schmidt since (10) is fulfilled 
for a > 1/2, compare Example 2.2.) Take an arbitrary orthonormal basis 
(ipk)keN °f the Cameron-Martin space of W\, and define f(v) and B(v) for 
v € Hq(0) analogously to (17). Then, Eq. (1) is an abstract version of prob¬ 
lem (16), and the mappings / and B are globally Lipschitz continuous (and 
thus linearly growing) from D((— A) 1 / 2 ) = Hq(0) to D((—A)°) = 1/2(0) 
and from D((— A) 1 / 2 ) to £(^ 2 ; D((— A) 1 / 2 ^)), respectively. The mapping 
properties of B follow from the fact that the Cameron-Martin space of W 1 
is continuously embedded in D((—A) l / 2 ~d) and the inequality 

W^Wd/z-a) 1 / 2 -P) ^ C\\ v \\ H i(o) IMId^-a) 1 / 2 -/ 3 )- 
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The latter is due to the inequalities ||utu||L 2 ( 0 ) 11 17 11 jyi (c>) 11 1(7 11 z , 2 (c?) together 
with 11vu;11< C'||u||#i( 0 )||iy||#i( 0 ) (a consequence of the Sobolev em¬ 
bedding H l (0) Loo^O) in dimension 1) and interpolation since 

D((-A)1/2-i5) = [L 2 (0),D((-A)^)\^ w . 

Thus, Assumptions 2.1, 2.3 and 2.5 are fulfilled for g = a = 1 / 2 , 0 < f3 < 1/4 
and 1/2 < a < 1 — 2 (3. 

Now let O C W l be d-dimensional and consider problem (16) where 
uq € Hq( 0), g : M —)• R is globally Lipschitz continuous, h : R —>• M is con¬ 
stant, and W\ = (lT r i(t)) tg [o i r] is a Wiener process whose Cameron-Martin 
space is continuously embedded in D ((—for some /3 < 1/2 — d/4. 
Then, the mappings / and B, defined analogously to the one dimensional 
case, are globally Lipschitz continuous (and thus linearly growing) from 
D^-A) 1 / 2 ) = Hq(O) to D({-A)°) = L 2 {0) and from D^-A) 1 / 2 ) to 
C(l 2 ] D{{— A) 1 / 2- ^)) respectively. It follows that Assumptions 2.1, 2.3 and 
2.5 are fulfilled for g = a = 1/2, f3 < 1/2 —d/4 and 1 < a < 1 — 2/3. As in Ex¬ 
ample 2.10 we could alternatively assume h : K —>• K to be sufficiently smooth 
and replace h(u(t,x )) in problem (16) by, e.g., h( f a k(x,y)u(t,y) dy) with 
a sufficiently smooth kernel k : O x O —> M. 

Remark 2.12. The reader familiar with SPDEs of the form (1) might won¬ 
der about the rather restrictive conditions in the examples above, especially 
on the noise terms therein. These restrictions are due to the fact that we 
basically adopt the setting from [35]. This allows us to focus on our main 
goal, i.e., the analysis of the error of the inexact counterpart of the Euler 
scheme (4), without spending too much time on explaining details regard¬ 
ing the underlying setting. However, it is worth mentioning that much more 
general equations of the type ( 1 ) have been considered in the literature, 
see, e.g., the recent results concerning the maximal L p -regularity of SPDEs 
in [33]. Also, the convergence of the linearly implicit Euler scheme has been 
considered under weaker assumptions, see, e.g., [12,13]. 

3 Exact Euler scheme 

In this section we use the linearly implicit Euler scheme to obtain a semidis¬ 
cretization of Eq. (1) in time. We present a corresponding convergence result 
as a slight modification of [35, Theorem 3.2]. Since no spatial discretization 
is involved, we speak of the exact Euler scheme in contrast to the inexact 
perturbed scheme considered in the forthcoming section. From now on, let 
Assumptions 2 . 1 , 2.3 and 2.5 be fulfilled. 

For K e N and t := T/K we consider discretizations (iqJjfL 0 given by 
the linearly implicit Euler scheme (4), i.e., 

u k+ 1 := (I - t A)~ l (- u k + rf(u k ) + y/rB(u k )xk ), k = 0,..., K — 1. 
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We use the abbreviation 


Xfc := Xk ■= ^ {W{t% + 1) - W{1 %)) 

with 

tk := r/c, k = 0,..., K. 

Note that each \ki k = 0,..., K — 1, is an -measurable Gaussian white 
noise on i 2 , i.e., a linear isometry from £2 to L 2 (Q, ^t k+1 , IP) such that for 
each a £<2 the real valued random variable Xk ( a ) is centered Gaussian with 
variance ||a||| . Moreover, for each k = 0,... ,K — 1, the sub-u-field of & 
generated by (xfc(a) : a £ i 2 } is independent of &t k - 

We explain in which way the scheme (4) has to be understood. Let G 
be a separable real Hilbert space such that D((—A) g ~P) is embedded into 
G via a Hilbert-Schmidt embedding. Then, for all k = 0,..., K — 1 and for 
all &t k -measurable, D((— A) e )-valued, square integrable random variables 
v € L2(S4, J^jP; D((— A) 0 )), the term B{y)xk can be interpreted as an 
-measurable, square integrable, G-valued random variable in the sense 

J 

B(v)xk := L 2 (n,^ t ¥]G)- lim ^ Xfc( b j) s («) b j, (18) 

J —^OO - 

J = 1 

where ( b j)j'GN is an orthonormal basis of £ 2 . This definition does not de¬ 
pend on the specific choice of the orthonormal basis Note that the 

stochastic independence of {xfc( a ) : a £ i 2 } and &t k is important at this 
point. We have 

nB(v)Xk\\h = nB(v)\\ 2 c Hs{ e 2 . G) , ( 19 ) 

the last term being finite due to the Lipschitz continuity of B by Assump¬ 
tion 2.3 (see also Remark 2.4) and the fact that the embedding 

D((-A) e ~P) G 

is Hilbert-Schmidt. Let us explicitly choose the space G in such a way that 
the terms u k + Tf(u k ) + y/rB{uk)xk on the right hand side of (4) can be 
considered as a G -valued random variable and the application of (/ — rA) _1 
to elements of G makes sense. Our choice of G, which we keep throughout 
this paper, is 

G := L)((-A)^ _max (°’ fT ’ /3+Q ' /2 )). (20) 

Note that the condition Tr(— A)~ a < 00 in Assumption 2.1 yields that the 
embedding D((—A) g ~P) D((—A) g ~^~ a / 2 ) is Hilbert-Schmidt, and the 

embedding D((—A) e ~P~ a t 2 ) D{{— J 4 )£’- max ( 0 >°'./ 3 +“/ 2 )) j s dearly continu¬ 
ous. Thus, we have indeed a Hilbert-Schmidt embedding D((—A) g ~P) G. 
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For all k = 0,..., K — 1 and v € L 2 (Ll, D((—A) g )) we consider the 

term B(v)xk as an element in the space 

L 2 (n,^ tk+ 1 ,¥-G) = L 2 (n,^ tk+ 1 ,¥;D((-Ay- m ^°^ +a ^)). 

Next, due to the Lipschitz continuity of / by Assumption 2.3 (see also 
Remark 2.4), we also know that for all v € L 2 (Q, LFt k , P; D((—A) e )) the 
term f(v) is an element in L 2 (Q, , P; G). Finally, as a consequence of 

Lemma 3.1 below and the fact that max(0,<r, /3 + a/2 ) < max(cr, 1/2) < 1 
due to (11), the operator 

(I-tA)- 1 : G D((—A) e ) 

is continuous. It follows that the discretizations (u k )\f =0 are uniquely deter¬ 
mined by (4) and for all k = 0,..., K we have 

Uk € L 2 (Q, JF tk ,F; D((-A) 0 )). 


Now we prove the lemma we just used to show the boundedness of the 
resolvents (I — rA) _1 of A in the right spaces. It will also be employed to 
prove Proposition 4.4 in the next section. 

Lemma 3.1. Let r > 0 and r € M. The operator I—tA is a homeomorphism 
from D((—A) r ) to D((—A) r ~ 1 ). For n €. N we have the following operator 
norm estimates for (I — TA)~ n , considered as an operator from D((—A) r ~ s ) 
to D({-A) r ), s < 1; 


\\(I - T A) "|| C(D((-A) r - S ),D((-A) r )) 



[(-Ar) s (l -rAi)-" 


0 < s < 1 
s < 0. 


Proof. The bijectivity of I — tA : D((—A) r ) —> D((—A) r ~ l ) is almost obvi¬ 
ous. Its proof is left to the reader. The bicontinuity follows from the conti¬ 
nuity of the inverse as shown below (case s = 1) and the bounded inverse 
theorem. Concerning the operator norm estimates, we use Parseval’s identity 
and the spectral properties of A to obtain 

sup _ \\(I - TA)- n v\\ 2 D{( _ A)r) 

ll t ’llr)((-A )’'- s ) =1 

= sup \\{I-TA)- n {-A) s - r wf m _ AY) 

\\w\\u=l 

= sup Wi-Ani-rAy^-Ay^wWl 

\\w\\u = l 

= sup ^2\(-X k ) s (l - T\ k )~ n (w,e k ) u \ 2 . 
fceN 
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If s < 0, the last expression is equal to (—Ai) 2s (1 — tAi) 2n . If 0 < s < 1, 
an upper bound is given by the square of 


sup 


*> 0(1 + TXY 


= S 6 ( 1-- 
n 


(n-s) 


(nr)' 


□ 


After these preparations we present an extension of the error estimate 
in [35], 

Proposition 3.2. Let Assumptions 2.1, 2.3 and 2.5 be fulfilled. Let (u k )k = o 
be the time discretization of the mild solution (it(t))t e [o,T] -D((— A) e ) to 

Eq. (1), given by the scheme (4). Then, for every 


6 < min(l — a, (1 — a )/2 — f3), 


we have for all 1 < k < K 

( E ||«(£fc) — w fe|||)((_ J 4)e)) — ^ ^ r<5 + ^(^II m o|Id((_A)£?)) 

where the constant C > 0 depends only on 5, A, B, f, a, /3, a and T. 

Proof. We argue as in the proof of Proposition 2.8 and consider the equation 

du(t) = ( Au(t ) + /(u(t))) dt + B(u(t)) d W (t), u(0) = uq, ( 21 ) 

where the operator A : D{A) C U —>• U is defined by U := D((—A) Q ), 
D{A) := D((— A) e+1 ), and Av := Av, v € D(A). Eq. (21) fits into the 
setting of [35], and its mild solution u = (u(£)) te ro ) Ti coincides with the mild 
solution to Eq. (1), compare the proof of Proposition 2.8. For K € N let 
(Ak)k =o b e gi ven by the linearly implicit Euler scheme 


Ho — UOi 

Uk+ 1 = (I - rAy 1 (u k + rf{u k ) + VrB{uk)xk) , k = 0,..., K - 1. 

By Theorem 3.2 in [35] we have, for all 5 < min(l — cr, (1 — a)/ 2 — /3), 

( E IH4) - u k \\l) 1 < c (V + ^(e||« 0 ||J) 7 ) , 

1 < k < K. The proof in [35] reveals that the constant C > 0 depends 
only on 5, A, B, /, a, /3, a and T. The assertion of Proposition 3.2 follows 
from the fact that the natural extensions and restrictions of the operators 
(I — rAy 1 and (I — rAy 1 to the spaces D{{— A) s ) = D{{— A) s+e ), set, 
coincide, so that u k = u k for all 0 < k < K, K G N. □ 
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( 22 ) 


Remark 3.3. (i) If k > K s (6 > 0), then l/k < T 5 t 5 , and we obtain 

(E||u(t fc ) — Wfe||x> (( - J 4 ) » ) ) ' < Ct& 

with a constant C > 0 that depends only on 6 , uq, A, B, /, a , /3, cr and T. 
Since 5 is always smaller than 1, it follows in particular that (22) holds for 
k = K , i.e., 

(E\\u(T)-u K \\ 2 Da _ A]e) y /2 <Ct 5 . 

(ii) The proof of Proposition 3.2 is based on an application of Theorem 3.2 in 
[35] to Eq. (15). The reader might have observed that therein the parameter 
s, which corresponds to our parameter a, is assumed to be positive. However, 
a closer look at the estimates in the proof of Theorem 3.2 in [35] reveals 
that the result can be extended to negative values of a and s, respectively. 
Alternatively, one can argue that if a < 0 then D((—A)~ a ) is continuously 
embedded into, say, D((— A) -1 / 2 ), so that Eq. (15) fits into the setting of [35] 
if / is considered as a mapping from U = D((—A)°) to Z?((— A) -1 / 2 ). We 
refer to [13] where the Euler scheme for stochastic evolution equations is 
considered in a more general setting than in [35]. 

4 Error control for the inexact scheme 

So far we have verified the existence and uniqueness of a mild solution to 
Eq. (1) as well as the convergence of the exactly evaluated linearly implicit 
Euler scheme (4) with rate 5 < min(l — a, (1 — a)/2 — /3). We now turn to 
the corresponding inexact scheme. We assume that we have at hand a ‘black 
box’ numerical scheme, which for any element w € L 2 (£l, &t k , IP; A) 0 ) 
approximates 

v := (I - tA)~ 1 (w + f{w ) + B(w)xk ) 

with a prescribed tolerance e > 0, the error being measured with respect to 
the L 2 (£l, ^t k+1 ! T; A) e )-norm. What we have in mind are w-wise ap¬ 
plications of some deterministic solver for elliptic equations, e.g., an optimal 
adaptive wavelet solver as developed in [9], combined with proper evalua¬ 
tions of the nonlinearities / and B, see, e.g., [11,16,30], and an adequate 
truncation of the noise B[w)xk■ We start with the initial condition uq and 
in each time step, we apply this ‘black box’ method to the approximation 
we have obtained in the step before. Our main goal is to prove that the 
tolerances in the different time steps can be chosen in such a way that the 
inexact scheme achieves the same convergence rate (in time direction) as 
its exact counterpart (Theorem 4.2). To this end we also analyze the error 
propagation of the inexact scheme (Proposition 4.3). 
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Our strategy relies on the ideas presented in [7]. Therein, Rothe’s method 
for deterministic parabolic equations is analyzed by putting it into an ab¬ 
stract framework, cf. in particular [7, Section 2]. We proceed similarly. There¬ 
fore, we first of all reformulate the exact linearly implicit Euler scheme (4) 
in the following way. We set 


n k := L 2 (Sl,& tk ,F-,D((-Ay)), 

o' 

II 

,,R, 

g k ■.= L 2 (n,^ tk ,r-G), 

k = 1, . 


Rr,k • T~Lk ^ Gk -\-1 

V R T ,fc(u) := v + r/(u) + y/r B(v)xk, 

o' 

II 

,,R- 

B T l '■ G k Rk 

v L t 1 v := (I — tA) l v , 

k = 1, . 

,,R. 


Recall that G = D((—A) e max (°> (7 ’P + a / 2 )) has been introduced in (20). With 
these definitions at hand, we can rewrite the scheme (4) as 

u k+1 := L^Rr^Uk), k = 0,... ,K - 1. (24) 

Remark 4.1. Without additional assumptions on B or a truncation of the 
noise expansion (18), the operator R r k cannot easily be traced back to a 
family of operators 


Rr,k,u) '■ D— A)^) —> G , uj £ fi, 

in the sense that for v £ R k = £ 2 (^ 2 , , P;-D((—A) e )) the image R r ,k(v ) 

is determined by 

(R t x(v))(uj) = R t ^,u}{v{u>)) for P-almost all uj £ fh (25) 

However, this is possible, for instance, if for all v £ D((—A) e ) the oper¬ 
ator B(v) : 1 2 —>• Z?((— A) e ~P) has a continuous extension B(v) : Uq —>• 
D((—A) e ~P) to a bigger Hilbert space Uq such that t 2 is embedded into Uq 
via a Hilbert-Schmidt embedding. Another instance where a representation 
of the form (25) is possible is the case where the mapping B : D((—A) e ) —>• 
C(£r,D((-A)^)) is constant, i.e., the case of additive noise. We take a 
closer look at the latter case, writing B € £(^ 2 ; D((— A) e_/3 )) for short. We 
fix a version of each of the P-almost surely determined, G-valued random 
variables Bxk = ^X k ■ & = 0,..., K — 1, K £ N, and set 

Rr,k,u{v) -=v + f(v) + (Bxk){u), uj € n, v £ D((-A) e ). 
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It is clear that (25) holds for all v € L 2 (Ll, D((—A) e )). Thus, in the 

case of additive noise, by setting 

n k :=D((-Ay)), 

g k -=G = D((-Ay~ max( - 0 ’ a ’ l3+a / 2) ), 

Rr,k,uj • R-k t Qk- 1-1 , . 

V H- R T ,k,u( v ) -= v + T f ( v ) + Vr(Bxk){u), 

L r 1 : Q k ~> Rk 

v i-)- Lf 1 v := (I — tA)~ 1 v, 

k = 0,..., K — 1, we can rewrite the scheme (4) in an cu-wise sense as 

Uk+i(u) := Lg l R r k)U (uk{u))■, k = 0,... ,K - 1. 

In this section we are focusing on the inexact counterpart of the scheme 
(4), which we introduce now. We assume that we have a numerical scheme 
which, for any k = 0,... , K — 1, any w € Rk-, an d any prescribed tolerance 
e > 0, provides an approximation [u] e of 


L r R Tj k(w), 


v = 


such that 



Given prescribed tolerances e k , k = 0,..., K — 1, for the spatial approxima¬ 
tion errors in each time step, we consider the inexact linearly implicit Euler 
scheme, defined as follows: 



(27) 


Note that the errors at each time step accumulate due to the iterated appli¬ 
cation of the numerical method [-] e . 

Next we present the main result of this paper. It describes a way how 
to choose the tolerances in the different time steps so that the inexact 
scheme (27) has the same convergence rate as its exact counterpart (24). 
Theorem 4.2. Let Assumptions 2.1, 2.3 and 2.5 be fulfilled, (u(t))t^[o,T] be 
the unique mild solution in D((—A) e ) to Eq. (1) and let 


5 < min(l — a, (1 — a)/2 — (3). 


If one chooses 
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for all k = 0,..., K — 1, K £N, then the output uk of the inexact linearly 
implicit Euler scheme (27) converges to u(T) with rate 5, i.e., we have 

1 /2 

(E\\u(T)-u K \\ 2 Da _ A)g) ) < Ct s 

with a constant C depending only on uq, 5, A, B, f, a, jd, o and T. 

Our strategy for proving this theorem relies on two auxiliary results, 
which we prove first. We start with Proposition 4.3, which states that the 
error propagation of the inexact linearly implicit Euler scheme can be de¬ 
scribed in terms of the Lipschitz constants C^f l J k of the operators 

Er,j,k = {L t 1 R T ,k- 1) ° (L r l Rr,k- 2 ) 0 ■ ■ ■ O (L r : Rj ~► Rk, (28) 

1 < j < k < K, K € N. Afterwards, we prove that these Lipschitz constants 
are bounded from above uniformly in 1 < j,k < K, K E N, see Proposi¬ 
tion 4.4. Finally, at the end of this section, we draw the proof of the main 
Theorem 4.2. 

Proposition 4.3. Let Assumptions 2.1, 2.3 and 2.5 be fulfilled. Let (u(f)) tg [ 0j T] 
be the unique mild solution in D((—A) e ) to Eq. (1). Let ( Uk)J? =0 be the dis¬ 
cretization of («(i))te[o,T] in time and space given by the inexact linearly 
implicit Euler scheme (27), where k = 0,..., K — 1, are prescribed toler¬ 
ances for the spatial approximation errors in each time step. Then, for every 
1 < k < K, K G N, and for every 

5 < min(l — a, (1 — a )/2 — jd) 

we have 

(e||u(4) - u k \\ 2 Da _ A)g) ^j 1 < C^T S + i(E||u 0 ||2, (( _ A)e) ) 

+ X C rJ+l,k £ j 
j =0 

with a constant C that depends only on 5, A, B, f, a, (3, a and T. 

Proof. We argue following the lines of the proof of [7, Theorem 2.16]. To 
this end, we introduce the operators 

Rr,k,k +1 • Rk t 

V 1 t E Ti k t k- Tl(^) •— [L r -R Ti fc(u)] £fc 
for k = 0,..., K — 1, and for 0 < j < k we set 

f-, J E r k—i,k o .. .o i , j k 

U j-u. 
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Note that, in order to define the operators E T j t k and E t j^ properly, we 
have to consider the measurability of the elements in their domains and 
in the corresponding image spaces, cf. our remarks concerning the right 
understanding of the scheme (4) at the beginning of Section 3. Therefore, in 
contrast to the situation in [7, Section 2], these Hilbert spaces depend on the 
indexes j and k. Nevertheless, with the right changes, we can argue along 
the lines of [7, Theorem 2.16]: We rewrite the difference Uk — Uk between 
the output of the exact and inexact scheme at time tk as an appropriate 
telescopic sum, so that by consecutive applications of the triangle inequality 
in Tik we obtain: 

I \u(t k ) - u k \\u k < | \u(t k ) - u k \\n k + II uk - uk\\n k 
< || u(tk) - Uk\\u k 

k- 1 

T ^ ^ ||^'rj,fc£'r,0,j('^o) E Tt j-\-i t kET,0,j+l(uo') ||^ • 

3=0 

Note that C^J k < oo for all l<j<k<K,KeN, because of Assump¬ 
tion 2.3 on the Lipschitz continuity of the free terms / and B and because 
of Lemma 3.1 on the boundedness of the resolvents of A. Thus, each term 
in the sum on the right hand side can be estimated as follows 

11 (^o) E T ^jj r \ kErfi,j+l M 11 ^ 

| ^T,j-\-l,kEr,j,j+l-E T fij (uo) E t j_\_i, k Er,0,j+l {uq ) 11 ^ 

— C'r/+l lfc ||^+l^r,0 J -(«0) — ^r,0,j+l( w o)||^. +1 - 

Since Uj = £ , Ti oj(uo), we obtain 

11 ^T,j,j+iE r ,Q,j (Ao)- E T fi j +1 («o)|| H . +1 

= \\Er,j,j+l{Uj) — E T jj + i(Uj)\\ n .^ < £j. 

Putting these estimates together yields: 


k -1 

IK4) - u k \\ Hk < || u(t k ) ~ u k \\ Hk + C'rJ+i ,k £ r 

3=0 

The error of the exact Euler scheme at time t k appearing on the right hand 
side can be estimated by using Proposition 3.2, which yields the assertion. 

□ 

In order to prove the main Theorem 4.2, it remains to verify the uni¬ 
form boundedness of the Lipschitz constants Y'Jk °f the operators E rj ^ 
introduced in (28). The proof is based on a Gronwall argument. 
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Proposition 4.4. Let Assumptions 2.1, 2.3 and 2.5 be fulfilled. There exists 
a finite constant C > 0, depending only on A, B, f, a, (3, a and T, such 
that the Lipschitz constants of the operators -E r ,j,fc introduced in (28) 

satisfy: 

<C for all 1 < j < k < K, I< e N. 

Proof. Fix 1 < j < k < K and observe that, by induction over k, 

Er,j, k(v) = L~^V 

k—j —1 

i =0 

for all v € TLj, where we set Ejj = I. Therefore, for all v, w G TLj , we have 


I E T .j } k ) E T j^{ui) \\'H k 

< \\Lf^v - Lf^wWn, 
k—j —1 

+ E r 


+ 


i=0 
fc-j-1 


U! 


H k 




i=0 


=:(/) +(IT)+ W- 


(29) 


We estimate each of the terms (/), (T) and ( M ) separately. 

By Lemma 3.1 and the trivial fact that ||u — w ||% fe = \\v — w\\y_, for all 
v, w € "Hj, we have 

(I) < H-C'r ll i C(D((-A)e))ll' l; ~ W \\Hk 

< (1 — tAi) - ^ - ^||u — w \\ u k (^0) 

< ||t> - W\\uy 

Concerning the term {E) in (29), let us first concentrate on the case 
a € (0,1). We use the Lipschitz condition on / in Assumption 2.3 and 
Lemma 3.1 to obtain 


w< 


k-j-l 

E T\\L^ k ~^ + \-AY 


1=0 


C(DK-Ap)) 


< 


< 


x Cj P || E Tt jj+i(v) E T jj+i(w )||^. +i 


k-j-l 

E - 


i =0 


MRll^C I—Ik 


k—j —1 

C/ L ' p E 


i=0 


/^iLip || || 

(r(fc — j — i)) a T dd+ l V W 


( 31 ) 
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For the case that a < 0 we get with similar arguments 

k-j-l 


(j) < cy p y 


r (~^l) <T /-(Lip 


i =0 


(1 — rAi) n r ’ j ' j+i 


(Ilu - w I 




< C / Llp (-A 1 )-’ £ TC“!A +i 

i=0 


( 32 ) 


v — w\ 


Let us now look at the term ( M ) in (29). Using the independence of the 
stochastic increments Xj+i and the equality in (19), we get 


k—j—i 


(M) 2 = Y te|| L^ k -^(B{E Tdd ^(v))-B{E TdJ ^(w)))xj+i 

i=0 
c-j-l 

X) r ll^ 


2=0 
k—j—1 

< -rll r-( fc -i)+dl 2 

i=0 


D«-A)e) 


C(D{{-AY-P- a /' 2 ),D((-A)e)) 


x E 


(-LrJJ+i^)) B ( y E T) jjj r i{w)') ^ Xj'+i 


D((-A)e-P- a / 2 ) 


k—j—1 

< Y, T \\ L T (k ~ j)+i \\ 2 c 

i=0 


£(D((-A)f?-/3-“/ 2 ),£>((-A)U) 


X E ||l?(£Vjj + j(u)) B(E T jj + i(w)) || J c Hg ( fe; r»((_ J 4)e-fl-o/a)) ■ 

Concentrating first on the case j3 + a/2 > 0, we continue by using the 
Lipschitz condition on B in Assumption 2.3 and Lemma 3.1 to obtain 

(J) 2 < *£ ’‘'r , J? + g/2) -(ll rv(-A)- 


i=0 


(r(fc — j — i)) 2 / 3+c 


x (d 1 ^ P ) ^ Il-®rj,i+i(' i; ) B T jj-\-i(w) II ZP((—v4)e) 


< (C]f) 2 Tr(-A)- 


(33) 


x 


fc-j-i 

X (r(fc - j - i)) 2 P+ a 






In the case (3 + a /2 < 0 the same arguments lead to 

(We 


i=0 


X (Cs P ) ® ll-^'Tj,j+i(' y ) -®'r,jj+*(' u; )||^((_A)e) 

<^ ip ) 2 Tr(-A)-“(-A 1 ) 2 ^“ £ r(C T L £ +i ) 2 ||u - 


( 34 ) 


t=0 
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Now we have to consider four different cases. 

Case 1. a £ (0,1) and /3 + a/2 € (0,1/2). The combination of (29), 
(30), (31) and (33) yields 


C Lip , < 

T,3,k — 


k—j — 1 


1 + C / iP E 


i =0 


T 

(' r{k - j - i)) a 


/''Tip 


+ C^ ip (Tr(-A)- 


«)l/2 


/k—j—1 

E 

V i=0 


1/2 


(r(ife - j - *)) 2/3+( 


(n Li P )2 


(35) 

Next, we estimate the two sums over i on the right hand side of (35) via 
Holder’s inequality. Set 


1 


q : = 


+ 2 > 2 . 


min(l — a, (1 — a)/2 — (3) 

Holder’s inequality with exponents q/(q — 1) and q yields 


k—j—i 

£ 

i=0 


c. 


Lip 


(r(fc - j - z)) ff T ’/d+i 


< 


/k—j—l 

E 


r 


0-1 

q /k—j—1 


q -1 


E 


i=0 


< 


T 


z' if 

E, 

Ki=l 

<1 l T t~ft it 


1 

k-j-1 \ o 

E 

i=0 / 


(36) 


L 


/k-j-i 


E 


Lip \<j 
j+i' 


i =0 


where the integral in the last line is finite since 2I2L = _ = 1. 
Similarly, applying Holder’s inequality with exponents (//(</ — 2) and q/2, 


fc-j-i 

E 

i =0 


(C 


Lip \2 


(r(k - j - i)) 2 ^+“ v T ’/b+i 


< 


< 


f k—j—l 

E 


_ (2/3 + a)q 

i=o (r(/c — j — *)) i -2 


<?— 2 2 

1 (k-j-l \ q 

E -(cffe w 

\ i =0 




(37) 


L 


q -2 

T (2/3+a)q \ q 

t I” 2 dt 


/fc-i-i 

n I J 




t=0 


The integral in the last line is finite since i—1. 


(2/3+a)q _ (2/3+q) (2/3+q) _ 

g-2 1-2/9 ^ 1—(l—o—2/3) 
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Combining (35), (36), (37) and using the equivalence of norms in M 3 , we 


obtain 


{C. 


T,j,k> 


<c 0 


k—j — 1 

1 + X t(G 

i =0 


Lip 



(38) 


with a constant Co that depends only on A, /, B, a , /5, a and T. Since (38) 
holds for arbitrary K E N and 1 < j < k < K, we can apply a discrete 
version of Gronwall’s lemma and obtain 


{CrJ k ) g < e (fc - j)TCo C 0 < e TCo C 0 . 

for all l<j<k<K,KeN and r = T/K. It follows that the assertion of 
the proposition holds in this first case with 

C := (e TC °C 0 ) 1/? . 


Case 2. (T < 0 and /3 + a/2 < 0. A combination of (29) with (30), (32), 
and (34) leads to 

C£& < 1 + Cj-'-f-A,)' £ r C^ + . 

i =0 

^_ j _| 

+ C^ ip (Tr(-7l)-“) 1 / 2 (-A 1 )^/ 2 ( T(C^ j+i ) 2 ) 1/2 . 

i =0 

Applying Holder’s inequality with exponent t /2 := 2 to estimate the first sum 
over * on the right hand side, we get 


k-j-l 


+ £ rfC^-u) 2 ) 1 ' 2 

i =0 

+ ^ ip (Tr(-4)- Q ) 1 / 2 (-A 1 )^/ 2 ( 2 r ^ L & +i ) 2 ) 1/2 > 


i =0 


which leads to 

(c“y 2 < c(i + £ T(cf& + ,) 2 ), 

t=0 

where the constant C E (0,oo) depends only on A, /, B, a, /3, cr and T. As 
in Case 1, an application of Gronwall’s lemma proves the assertion in this 
second case. 

Case 3. a € (0,1) and /3 + a/2 < 0. In this situation, we combine (29) 
with (30), (31) and (34) to get 


^Lip 

r,j,k 


< 


k-j-l 


1 + c / ip E 


i =0 


T 

(r(fc - j - i)) CT 


/^Lip 

j'+i 


fe-j-1 


+ C^ p (Tr(-A)- Q ) 1 / 2 (-A 1 ) /3+a/2 ( X] 


»=o 


1/2 


24 



Setting 


1 


Q3 ■ = 


1 — a 


+ 2 


and following the line of argumentation from the first case with q% instead 
of q we reach our goal also in this situation. 

Case 4. a < 0 and (3 + ot /2 E (0,1/2). Combine (29), (30), (32) and (33) 
to get 


c% t < i+ cy-i-Mr E tC tZ-+ i 

i=o 

+ C L B ip (Tr(-A)~ a ) 1/2 ( 

V i=0 


1/2 


(r(fe -j- i)) 2 P+ a v T " +i 


(G 


Lip \2 


Arguing as in the third case with 

1 

Q4 '- 1/2-((5 +a/2) + 2 

instead of q 3 , we get the estimate we need to finish the proof. □ 

We conclude with the proof of our main result. 

Proof of Theorem f.2. The assertion follows from Proposition 4.3 combined 
with Proposition 4.4 by using the elementary estimates 

1 1 K ~ 1 

— < —j = T~ s t s and Ve t < Tt 5 . □ 

K - K 5 ^ ~ 

j =0 
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